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We consider a two-dimensional electron gas in an inversion asymmetric layer and in the presence of 
spatially distributed magnetic impurities. We investigate the relationship between the geometrical 
properties of the wave-function and the system's spin dependent transport properties. A localiza- 
tion transition, arising when disorder is increased, is exhibited by the appearance of a fractal state 
with finite inverse participation ratio. Below the transition, interference effects modify the carrier's 
diffusion, as revealed by the dependence on the scattering time of the power law exponents charac- 
terizing the spreading of a wave packet. Above the transition, in the strong disorder regime, we find 
that the states are spin polarized and localized around the impurities. A significant enhancement 
of the spin current develops in this regime. 

PACS numbers: 72.25.Dc, 72.25.Rb, 75.50.Pp 



I. INTRODUCTION 

The intrinsic spin Hall effect refers to the generation 
of a dissipationless spin current in response to an ex- 
ternal transverse electric field, it was predicted to occur 
in p-type semiconductors 1 and n-type heterostructures, 2 
and it was observed experimentally in non-magnetic two- 
dimensional systems. 3-5 The effect originates in the cor- 
relation between the particle's motion, driven by the ex- 
ternal electric field, and its spin, through the spin-orbit 
interaction. In the case of a two-dimensional electron gas 
with spin-orbit Rashba coupling (in an inversion asym- 
metric layer), the spin Hall conductivity is suppressed by 
impurity scattering. 6 However, breaking time invariance 
by the presence of magnetic impurities, restores the spin 
Hall effect. 7 9 Remarkably, spin transport in a magnetic 
disordered layer is reinforced by multiple scattering inter- 
actions, and, for a large range of impurity concentrations 
(in the weak disorder regime), the spin conductivity re- 
mains near its universal clean value, ct s h = — e/8ir (— e 
is the electron charge). 9 Two questions arise naturally 
concerning the behavior of the spin Hall conductivity: 
first, how does it evolve in the strong disorder regime, in 
particular, does the system undergo a localization transi- 
tion; second, what is the influence of the carrier mediated 
interactions between impurities, that may induce a ferro- 
magnetic transition. In this paper we shall focus on the 
first problem, the dependence of ct s h on the concentra- 
tion and strength of the paramagnetic disorder. 

The understanding of spin transport mechanisms in 
semiconductors, in which the existence of a strong spin- 
orbit coupling enables the electrical manipulation of the 
spin, is important for spintronics applications. 10 14 In 
particular, recent experiments in magnetically doped 
semiconductors, 15 ' 16 reveal interesting physical phenom- 
ena, such as the generation of spin currents, 17 the fractal 
geometry of localized states, 18 and intriguing properties 
of the intrinsic anomalous Hall effect. 19 ' 20 These exper- 
iments show the interplay between confinement, spin- 
orbit splitting, magnetic disorder and carrier mediated 
interactions. Their combined action drives spin trans- 



port, determines the magnetization structure, and gov- 
erns the localization transition, that modifies in turn the 
magnetic properties and the spin conductivity. A mini- 
mal model Hamiltonian accounting for these effects, must 
contain a hopping term, a Rashba term coupling the car- 
rier's momentum to their spin, and an exchange term 
with the interaction of the itinerant spins with the mag- 
netic moments of randomly distributed impurities. 21 24 

Most studies on dilute magnetic semiconductor quan- 
tum wells, deal with the influence of the carrier induced 
ferromagnetic transition and phase separation on the 
transport mechanisms, and more specifically, with the 
existence of carrier localized states and their impact on 
the anomalous Hall effect. 25 As a natural extension to 
these investigations, we focus in this paper on the in- 
trinsic spin Hall current, which arises, as the analogous 
intrinsic anomalous Hall conductivity, from the Berry 
phase contribution to the velocity operator. 26 Our goal 
is to identify the Anderson localization of the carriers 
as a function of the disorder strength, and to relate the 
geometry of the quantum states to the spin transport 
properties beyond the weak disorder regime. 

We start with a one band tight binding model taking 
into account spin-orbit coupling and impurities interac- 
tions. In this framework we study the spreading of a 
wave packet by solving the time-dependent Schrodinger 
equation. The knowledge of the time evolution allows the 
computation of the typical diffusion exponents, the frac- 
tal correlation dimension and the inverse participation 
ratio that characterize, at long time scales, the transport 
and localization of the quantum states. We find that 
the exponents vary from near ballistic values in the weak 
disorder regime, to a diffusion limit at stronger disorder. 
Increasing the disorder, one observes that the wave func- 
tion becomes fractal and that the inverse participation 
ratio remains finite, indicating the existence of a local- 
ization transition. This is further confirmed by the local 
density of states that concentrates around the impurities. 

The presence of interference effects in the spreading 
of the wave function are investigated, in the weak disor- 
der approximation, using the linear response theory. We 
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compute the quantum corrections to the spin conduc- 
tivity, and find that a weak antilocalization effect arises, 
which is essentially proportional to the spin-orbit induced 
spin splitting. 

The spin Hall conductivity, which initially decreases 
with the disorder strength, undergoes a drastic increase 
in the strong disorder limit, especially when impurity 
bands appear in the density of states. This enhance- 
ment of the spin conductivity in a regime where the car- 
rier states are localized is attributed to the tight correla- 
tion between the carrier and impurity spins, that locally 
breaks the symmetry between spin up and down states. 

II. MODEL 

We consider an electron conduction band of a semicon- 
ductor heterostructure, doped with a concentration c of 
magnetic impurities, in the tight binding approximation. 
We assume a simple geometry, where electrons and im- 
purities reside in a square lattice of TV sites, size L x L 
and spacing a. The model hamiltonian, 

H = H e + H so + H xc (1) 

contains a kinetic, 

a spin-orbit, 

H so = — ^2 (^l^Ci^y - c\a y Ci+£j + h.c. 

i 

and a double exchange term, 

H xc = -^yZ ni ' c l ac ii 

where a is the vector of Pauli matrices, c\ = (c\^ c\^) 
(ci) is the creation (annihilation) operator on site i and 
spin up t, or down ^; thop is the hopping energy between 
neighboring sites z, j, A is the Rashba spin-orbit coupling 
constant, and J s the exchange interaction energy con- 
stant. The spin-orbit term, proportional to A, contains 
hops combined with a spin-flip to neighbors in the x and 
y directions. The magnetic moments are randomly dis- 
tributed over sites i belonging to the set X of impurities, 
and their orientation n^, is uniformly distributed on the 
unit sphere. We choose units such that h = thop = a = 1. 

The above Hamiltonian, can be considered as a min- 
imal model of a dilute magnetic semiconductor electron 
gas, confined in a well between two asymmetric layers, in 
the paramagnetic state. 27 It is well suited for investigat- 
ing the anomalous Hall effect, and the intrinsic spin Hall 
effect in the presence of magnetic impurities. 28 

The presence of spin-orbit coupling and magnetic im- 
purities breaks spin rotation and time reversal symme- 
tries, allowing the occurrence in two dimensions, of the 



Anderson metal-insulator transition. 29 31 To characterize 
the localization of the quantum states and its influence 
on spin transport properties, it is convenient to inves- 
tigate both dynamical and spectral properties. 32 34 We 
therefore investigated the spreading of a wave packet and 
the spatial distribution of the density of states. From the 
time evolution of the wave function, we can extract in- 
formation on the diffusion and the time autocorrelation 
exponents, which are related with the spectrum of fractal 
dimensions. 35 ' 36 To distinguish between extended and lo- 
calized states we use the inverse participation ratio. 37 39 
The localization transition is usually associated with a 
mult ifr act al structure of the wave functions, that can 
deeply influence the transport mechanisms. 40 42 We also 
measure the local density of states to characterize the 
spatial distribution of the eigenstates over the energy 
spectrum. 43 ' 44 

The quantum state of the itinerant spins \ijj(t)), evolves 
according to the Schrodinger equation 

i^- t \m) = H\m), (2) 

which gives, in position representation, the spinor wave 
function ip(x,i) = (x\ijj(t)), where x belongs to the 
square lattice, and which spin up and down components 
are (^, ip±). We use the kernel polynomial method to in- 
tegrate (2). 45 ' 46 The initial condition is a Gaussian wave 
packet of size d and unit norm, 

^0) = ^i_exp (3) 

which spin is in the up direction ip = (^,0). Most nu- 
merical computations were performed using L = 2048 
and d = 8. Figure 1 illustrates the evolution of the wave 
packet, for two values of the disorder strength. 
The width of the wave packet is defined by, 

w 2 n (t)= [ dx\x\ 2 \Mx,t)\ 2 . (4) 

JL 2 

The asymptotic time dependence of the width, indepen- 
dent of the spin, satisfies a power law, 

w(t) - t p (5) 

characterized by the exponent /?; typical values are /3 = 1, 
for a free particle (the Rashba term do not change this 
linear dependence on time), and j3 = 1/2, in the diffusion 
regime (see Fig. 2). 

The geometry of the time dependent quantum state 
can be described by a measure based on the electron wave 
function ^(t), which gives the probability to find the 
particle at site z, |^(t)] 2 . Dividing the lattice into Ni 
cells of size Z, the dynamical measure of the k box is then 
given by, 

m k (l,t)= J2 I^WI 2 ' fc = l,...,JV z , (6) 

ieB k (l) 



FIG. 1. (Color online) Time evolution of the wave function; we use a logarithmic scale log \ip^(x, t)\, to visualize the small 
amplitude region. The lattice size is 2048 2 . The exchange constant is J s — 0.1 (top row), and J s — 1.0 (bottom row); times are, 
from left to right, 160, 288, 416, 544, 672, 800. The concentration of magnetic impurities is c = 0.1, and the Rashba constant 
A = 0.2. The color scale is fixed, and spans the interval (10 -4 , 5 10 -2 ). 



where we sum over the I 2 sites of the Bk box. Using this 
measure we define the dynamical inverse participation 
ratios, 



P q (l,t) 



k=l 



and the spectrum of fractal dimensions, 

l logP g (M) 



D q (t) lim 



i^o q - 1 log I ' 



(7) 



(8) 



which determine, in the long-time limit, the asymptotic 
behavior, 39 ' 47 ' 48 

P q (l) ~ lte-V D * . 

In the following we show results with q = 2, of the dy- 
namical inverse participation ratio and correlation di- 
mension, P2(t) and D2(t) (see Figs. 3 and 4). 

The local, spin dependent, density of states is defined 
by, 

p t (x,e) = Y,cos 2 (e/2)\^(x,a)\ 2 S(e - e a ) (9a) 

a 

Pl (x, e)=J2 sin 2 (0/2)|^(aj, a)\ 2 6(e - e«) (9b) 

OL 

for the spin- up and spin-down components, respectively, 
where is the angle of the impurity magnetic moment, 
chosen as the quantization axis, with respect to the z di- 
rection, e is the energy and e a is the eigenenergy of the 
eigenstate \a) (we use the notation i(j s (i,a) = (i\s,a), 
with s =t,l is the spin index). The sum over the space 
gives the total density of states p(e) = J]JPt(^ e ) + 
p±(i,e)]. As for the time dependent Schrodinger equa- 
tion, we use the polynomial kernel method to compute 
(9). 9 ' 46 This quantity is experimentally accessible, and 



shows a multifractal structure in the vicinity of the metal- 
insulator transition, on the surface states of a dilute mag- 
netic semiconductor. 18 The explicit dependence on the 
spin is useful for investigating the correlation of possi- 
bly localized electron spin states, with the spin of the 
impurities (see Fig. 5). 



III. WAVE PACKET DYNAMICS AND 
LOCALIZATION 

In the usual Anderson model one considers a free par- 
ticle that jumps between sites having random energies, 
independently of its spin; in the present model, a particle 
jumps changing its spin, due to the spin-orbit coupling, 
and is scattered off by spatially distributed impurities. 
The study of the dynamics of an initially localized wave 
packet, is interesting for getting a qualitative image of 
the propagation of the quantum states in the field of ran- 
domly distributed scatters. We show in Fig. 1, the time 
evolution of (3) for two values of the disorder, J s = 0.1 
in the top panel, and J s = 1 in the bottom panel. Most 
numerical results are given for a concentration c = 0.1 
of impurities, and a Rashba constant A = 0.2. Not only 
the speed of propagation differs in both cases, but also 
the geometry of the wave function is qualitatively differ- 
ent. In the weak disorder case, one observes an interfer- 
ence pattern, and at long times, a star-like structure; for 
stronger disorder, the quantum probability density de- 
velops a granular and intermittent structure, with a slow 
decay of the amplitude at the origin. 

These qualitative observations are confirmed more 
quantitatively by the measure of the wave packet's mean 
square displacement w(i), that characterizes its diffusion, 
and the correlation dimension, which is related to the 
decay exponent of the return probability. As shown in 
Fig. 2, the asymptotic behavior of the width is well de- 
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FIG. 2. (Color online) Width of the wave function showing 
a time power law with exponents dependoing on the disorder 
strength, from quasi free dispersion ()3 = 0.95, for J s — 0.1) to 
diffusion (/? = 0.51 for J s — 1.0). Exchange constant J s — 0.1, 
0.2, 0.4, 1.0, and 2.0 (from top to bottom). 
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FIG. 3. (Color online) Measure of D2 as a function of time 
for different disorder strengths (J s in the legend). 



scribed by a power law (5), with an exponent j3 that de- 
creases with increasing disorder. We find that for J s C 1 
the motion is quasi ballistic, while for J s > 1 a slightly 
subdiffusive regime sets in, as a possible manifestation of 
localization. The strong dependence of j3 on the disorder 
strength is reminiscent of the diffusion in a quasi periodic 
potential. 33 ' 49 ' 50 

The measure of the correlation dimension D 2l which 
accounts for the number of sites visited at time t, w(t) D2 , 
also depends on the disorder strength. In Fig. 3 we plot 
D 2 (t) for J s = 0.1, 0.2, 0.4, 1.0, 2.0, and 4.0; the corre- 
sponding asymptotic values of the fractal dimension are 



FIG. 4. (Color online) Inverse participation ratio P2 as a 
function of time for different disorder strengths (J s in the 
legend). When J s > 1, P2(t) remains finite for long times, 
indicating the existence of a localization state. 



D 2 = 1.45, 1.32, 1.14, 0.99, 0.88, and 0.79. It is worth 
noting that for J s = 1 we find /3 « 0.5 and D 2 « 1, these 
are just the exponents of a two-dimensional random walk, 
giving a strong indication of a critical state separating 
extended and localized regimes. For stronger disorder 
the correlation dimension passes below 1, D 2 < 1, corre- 
sponding to a point-like support of the wave packet, and 
a persistence of the probability to return to the origin, in 
accordance with the phenomenology observed in Fig. 1. 

A criterium to distinguish between extended and lo- 
calized states is given by the asymptotic behavior of the 
dynamical inverse participation ratio P 2 (t), that must 
vanish for extended states but remains finite for local- 
ized ones. We plot P 2 (t) in Fig. 4, using the same pa- 
rameters as before, and remark that a visible change op- 
erates around J s = 1; the asymptote is not only finite 
for J S >1, but the temporal fluctuations are notably en- 
hanced around the transition region. The picture shows 
a portion of the temporal evolution; for J s = 2 and 4, 
the simulations run up to 5 x 10 3 and 10 4 , respectively, 
without change in the mean value of P 2 . Therefore, from 
the dynamical point of view, a localization transition ap- 
pears when the disorder is increased, around the values 
c = 0.1 and J s = 1. 

In order to complete the dynamical picture of the quan- 
tum system with its spectral properties, we investigate 
the local density of state (9). We find that for weak 
disorder the (spin independent) local density of states 
p(i) = pi(i) + PtW? an d the total one p, coincide: there 
is no trace of strong inhomogeneities, for a given random 
position of the impurities. However, spin fluctuations are 
important even at weak disorder. Figure 5 presents the 
spin dependent local density of states, for the weak dis- 
order case, upper panel J s = 0.1, and for the transition 
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FIG. 5. (Color online) Spin up and down local density of 
states for J s = 0.1 (top row) and J s = 1 (bottom row). 
The contours (blue and red) are at the location of the im- 
purities, where the z component of their magnetic moment is 
rtz — —0.5 (blue) and —0.5 (red). In the neighborhood of the 
impurities, at weak disorder strength, ptW ano ^ com_ 
pensate to give a value close to the one of the total density 
of states p; at variance, for a stronger disorder, the spin of 
the quantum states is strongly correlated with the orienta- 
tion of the impurity magnetic moment. The Fermi energy is 
around 2, and the bottom color bar gives the scale of p\ only 
a fraction 32 x 32, of the total lattice is shown. 



region case, bottom panel J s = 1. When the disorder 
is weak, in the delocalized state, the spin fluctuations at 
the impurities sites compensate each other. Note in par- 
ticular, the exchange of colors between the up-spin and 
down-spin images, around the impurities contours (top 
row). In contrast, for stronger disorder (bottom row), a 
depletion of the density of states at the spin-up locations 
is created, while at the spin-down locations, an enhance- 
ment is observed. In the incipient localized state, for 
parameters near to the transition, arises a strong corre- 
lation between the spin of the electrons and the direction 
of the magnetic moment of the impurities. This effect 
on the spin translates in a strong inhomogeneity of the 
spin independent p(i): total and local density of states 
no more coincide, and as in the usual Anderson transi- 
tion, the arithmetic and geometric means of p(i) are no 
longer equal. 



IV. SPIN HALL CONDUCTIVITY 

At weak disorder the spin Hall conductivity stays close 
to its clean value, 9 but already quantum interference ef- 
fects play an important role (Fig. 1). This observation 
naturally leads to questions about the effect of interfer- 
ence on the spin Hall conductivity. Quantum correc- 
tions for a system with spin-orbit coupling and spin- 
independent scattering, were previously calculated for 
the charge conductivity in connexion with the anomalous 
Hall effect, 51 as well as for the spin Hall conductivity. 52 
Another question that arises from the investigation of the 
localization transition is that of the persistence of the 
spin Hall effect in the localized regime. The weak local- 
ization corrections, due to quantum coherent backscat- 
tering, are therefore of interest because they may give 
information about the quantum interference mechanisms 
that should also be relevant for spin transport in the 
strong disorder regime. 



A. Quantum corrections 

In order to understand quantum effects of the dis- 
ordered magnetic impurities on the spin Hall effect 
we will calculate the first quantum corrections to the 
Kubo formula for the spin Hall conductivity using the 
parabolic continuous band approximation of the tight- 
binding Hamiltonian (1). The continuous Hamiltonian is 
written 



H 



EL 

2m 



A 



(z x p) + V{x) , 



(10) 



where p and x are the momentum and position opera- 
tors of the electron in the plane (x,y). The interaction 
potential is characterized as in the tight-binding model, 
by the exchange energy J s and a characteristic micro- 
scopic length scale a (the lattice step in the tight-binding 
model). It can be written, 



V(x) 



COS Vi 



ri sin t 
cos Oj 



6(x — Xi) . 



(ii) 

where and Oi and <pi are the angles of the magnetic mo- 
ment of an impurity located at X{. The concentration 
of impurities will be denoted by c. We work in nondi- 
mensional units, such that h = rn = a = e = 1. The 
spin splitting at the Fermi energy is given by A = 2\kp 1 
where kp = is the Fermi wavenumber. The clean 

Hamiltonian is diagonal in the chiral basis, with eigenen- 
ergies 



e±(fc) 



k 2 



± Xk . 



(12) 



where k is the modulus of the wavevector k = (k x , k y ) — 
k (sin ip, cos ip). 

The linear (spin Hall) response of the system, to an 
applied electric field E X1 is given by the Kubo formula 
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FIG. 6. The three Hikami boxes for the weal localization 
corrections to the spin conductivity. The spin and charge 
current vertex are renomalized. 
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FIG. 7. Diagrammatic representation of the Cooperon equa- 
tion. 



for the spin Hall conductivity. It is convenient to express 
the Kubo formula in terms of retarded G R , and advanced 
G A , Green functions, as 52,53 



o* xy {u>) =- (Tr [f y (k)G R (e + aj,k)j x (k)G<(e,k) 

+ f y (k)G<(e + io,k)j x (k)G A (e,k)]), (13) 

where 

72 I 



™ r f d 2 k f de m 



(-) 



Tr s is the trace over the spin index, (...) denotes the 
average over the spatial disorder and magnetic moment 
orientations, and G < (e) = f(e)(G R (e) — G A (e)), where / 
is the Fermi function. 

The charge current operator is j x = — v x = — i [H, x] 
and the spin current operator is j* = (1/4) {cr^, v y }. The 
disorder averaged retarded and advanced Green functions 
are written 51 

where the spin-flip time r is defined by the self-energy as 
1 



2itN f cJ z s « cJ z s . 



We assume that the spin splitting is small and disor- 
der weak to that the Fermi energy is the largest energy 



scale of the system. This leads to the approximation 
€f ^> A,r -1 and Np w 1/(2tt). The nondimensional 
parameter x = At contains information about the spin- 
splitting and the disorder strength, and can take arbi- 
trary positive values. Using this parameter, the Green 
functions can be written in a nondimensional expression 
as 

nR,A ft tn \ _ -£(k)±i + x(a x suk/9 - a y cosy?) 

° { ^ ~ (-Z(k)-x±+i)(-Z(k)+x±i) ' (16) 

where £ = 2r(/c 2 /2 — €f)- The fundamental interaction 
vertex Tq is written in terms of the impurity concentra- 
tion and the interaction term as 



= c f V(0, 



(2ttN f )3t 



^)^(M): 



no o o^ 

0-120 
2-10 
\0 ly 



(17) 



where the four center elements are non-zero due to spin 
mixing. 

For convenience we write all expressions in matrix no- 
tation, using the following definitions: 

is the definition of the Kronecker tensor product and 
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is the definition of matrix multiplication. For the sake 
of clarity we will note 4x4 matrices using a calligraphic 
letter style. 

The spin Hall conductivity in the linear response 
framework, can be computed as a perturbation series in 
the disorder strength, 



+ 0" 



sH 



(18) 



where the first order takes into account the zero- 
interaction-loop spin conductivity (equivalent to the 
Drude term for the electrical conductivity), and the other 
ones contain the classical ladder corrections, plus lowest 
order quantum corrections. The first two terms in (18 
were computed in Ref. 9; here we compute the third term 
(maximally crossed diagrams). 

The terms in the perturbation series with the quantum 
corrections 5<j s h, can be diagrammatically represented 
by the so-called Hikami boxes, as shown in Fig. 6. 54 ' 55 
The total first order quantum corrections are the sum 



(15) of the contributions of the three boxes Sa s n = 8&sh 



Using the diagrams of Fig. 6, we obtain the following 
expressions for the three Hikami boxes: 



7 



S^H = ~ It (G A {k)J*G R (k)) ® (G R (q - k)J x G A {q - k)) C(q) , (19) 

J q J k 

^?h = -[ flf (G A (k)J*(k)G R (k)V(0,<j))G R (k')) 

JqJaJkJk' 

® (G R (q - k', )J x (q - k')G A (q - k')V(9, <P)G A (q - k))C{q) , (20) 
^?h = -[ fit (G A (k)V(0,cj ) )G A (k')J^k')G R (k')) 

J q JQ Jk Jk' 

® (G R (q - k')V(0, 4>)G R (q - k)J x (q - k)G A (q - k))C{q) , (21) 



where we used the notations J = J dq/ (27r) 2 , f n = 

f dO d(j) sin / '4tt is the integration over all magnetic mo- 
ments configurations. In (19)-(21), J x is the renormalized 
charge current vertex 



J X 



-A (2 + x 2 ) 



(22) 



where we remark the appearance of an extra A factor, 
and Jy is the renormalized spin current vertex of z-spins 
in the ^/-direction 



J* 



kx 



2(8 + 7x 2 ) 



(23) 



Taking into account the renormalization of the vertices 
turns out to be important, as they will rescale the quan- 
tum corrections to the spin conductivity with a factor 
proportional to A/ep. This correction originates in the 
anomalous contribution to the velocity operator due to 
the spin-orbit coupling. 

C is the Cooperon, a 4 x 4 matrix coupling the incom- 
ing to the outgoing spin states. Physically it gives the 
amplitude of corrections due to quantum interference of 
counter propagating trajectories. The Cooperon satisfies 
the self-consistent Bethe-Salpeter equation, that is most 
easily expressed in terms of the diagram in Fig. 7. In or- 
der to derive the spin conductivity we need the Cooperon 
at zero frequency (static) and at the Fermi energy; in this 
approximation it satisfies the following equation, 

C(q)=r +r I G R (p)®G A (-p + q)C(q). (24) 

The derivation of the explicit expression of the Cooperon 
is presented in the Appendix. The qr-dependence of its 
elements can be schematically expressed in the form [see 
Eq. (Al)] 



a(x,e^) (b(x)+l 2 q 2 ) 
2ttN f t (c(x) + l 2 q 2 )(d(x) + l 2 q 2 ) ' 



(25) 



where a, 6, c and d are nondimensional rational expres- 
sions of x and e 1% , and x is the angle between k and q; 



I = rkp is the mean free path. From the explicit expres- 
sion (Al) we obtain, the behavior of the Cooperon in the 
large x limit, 



lim C{q) = -±—: 



I 16(£V + 10) 
' 3Z 4 g 4 +80JV+400 




^i X] 2 q 2 





jV±2_ 



/-(/■■ l 2 q 2 -10 
6(7V-4) 6(7V-4) 
Z V-10 l 2 q 2 +2 



6(7 2 g 2 -4) 6(7V-4) 
\3^g 4 +80^q2 +4 oo u u 



-8e~ 2 ^V \ 
3Z 4 q 4 +80Z 2 q 2 +400 




16(^ 2 g 2 + 10) f 
3Z 4 q 4 +80£ 2 g 2 +400/ 

(26) 



For small values of x we get, 



lim C(q) = - Ar : 
x^o 2ttN f t 



( 4+l 2 q 2 




V o 









2l 2 q 2 


8 


-16-MV 


16-l 2 q 2 


8 


2l 2 q 2 


16-l 2 q 2 


-16-MV 












2 

A+l 2 q 2 



(27) 



From the above two expressions we note that the q —> 
limit of the Cooperon, changes according to the value of 
x. This implies that we must retain the full dependence 
on x in order to obtain a reliable result for the spin con- 
ductivity. Indeed, in the calculation of the spin Hall con- 
ductivity we are interested in the q close to zero limit, 
that allows us to simplify the qr-dependence of the Green 
functions and to factor out, in (19)- (21), the q integra- 
tion of the Cooperon. Because the interaction is isotropic 
on average (paramagnetic state), the two elements con- 
taining the angle x will vanish after integration. One can 
also see from these limits that to first order the Cooperon 
will give logarithmic terms in q. 

Details on the computation of Eqs. (19)-(21), are given 
in the Appendix. We find that the quantum corrections 
contribute positively to the spin Hall conductivity (effect 
usually referred to as weak anti- localization). The correc- 
tions of all three boxes are proportional to the spin-orbit 
splitting energy A/ep- We present in Fig. 8 the plot of 
Eqs. (A6) and (A7), that give the Hikami boxes correc- 
tions to the spin conductivity. The second and third 
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FIG. 8. (color online) Quantum corrections to the spin Hall 
conductivity, for the first, second and third (identical) Hikami 
boxes, and the total corrections. In units of (— e/87r)(A/ei?) 
the quantum corrections are function of x = At only. 



boxes, that contribute equally to the corrections, are 
about one order of magnitude smaller than the correc- 
tions of the first box, and their contribution is negative. 
For large values of the parameter [see Eqs. (A6-A7)], 



Kh = -- 

£ ( 2 ) 

sH = ~ 



1 A (81og(f)-71og(f)) 
87r e F 784 x 

1 A log (jf) 
8tt e F 3528 a; ' 



(28) 
(29) 



the corrections tend to zero as 
for small values of x, 



1/x. They also vanish 



-1 

8tt e^64' 



8aW 



-1 A f -x 
8tt V 192 



(30) 



In the intermediate region, for x of order one, the quan- 
tum effects reach a maximum. It is interesting to note 
that the quantum corrections, up to a global factor A/ep, 
are determined by a universal function, displayed in 
Fig. 8, that depends only on the ratio of the spin-orbit 
splitting to the disorder strength. The shape of this func- 
tion, shows that increasing the disorder (with fixed spin- 
splitting), the quantum contribution to the conductivity 
first increase, and then, after reaching a maximum, de- 
crease. 

It is worth noting that we calculate a spin conductiv- 
ity and therefore the term (anti-) localization does not 
exactly have the same meaning as in the charge conduc- 
tivity case. (Anti-) localization must be taken in the 
quantum spin transport sense and its meaning is more 
subtle than in the quantum charge transport picture. 



B. Spin conductivity 

Disorder has a very strong effect on the spin conduc- 
tivity, especially for high values of the exchange energy, 
when separate impurity bands are present. For weak dis- 
order we observe that the static spin conductivity slightly 
decreases with respect to the clean limit, 9 including a 
small modification due to quantum interference effects. 
However, at strong disorder the behavior of ct s h qual- 
itatively changes. In Fig. 9 we plot ct s h as a function 
of the Fermi energy, for different values of the disorder 
strength, from the weak to the strong disorder limits. In- 
creasing the disorder strength, for values of J s < 1, cj s h 
decreases, and for J s > 1, it increases, showing that the 
localization transition at J s « 1, modifies the behavior of 
the spin conductivity; this can be naturally attributed to 
the strong correlations of the carrier spin with the impu- 
rities, as discussed in Sec. Ill, resulting from the fractal 
geometry of the quantum states and their localization. 

If disorder is strong enough, the total density of states 
develops impurity bands; in this case, a striking intermit- 
tency in space and energy of the local density of states 
arises (Fig. 10). We observe that, near the localiza- 
tion transition, the energy distribution of states differs 
between sites, depending strongly on whether they are 
occupied or not by an impurity. For a given impurity, 
electronic states display a large distribution of energies. 
However, we observe important variations from site to 
site, showing a pronounced bias between positive and 
negative energy states. This systematic asymmetry is re- 
lated to the spin dependence of the localized states. For 
strong disorder, there are sites for which the quantum 
states have energies almost entirely concentrated in the 
impurity bands, these states are strongly spin polarized. 

The striking enhancement of the spin conductivity, 
starting at the localization transition, can be explained in 
terms of the strong spatial spin fluctuations at the impu- 
rity sites. Indeed, the spin current results from the imbal- 
ance between counter propagating spin-up and spin-down 
electrons; the presence of a magnetic impurity locally 
produces a marked bias in the spin orientation, due to 
the electron wave function localization: only states with 
the appropriate spin polarization are allowed. As we have 
shown before, quantum corrections, proportional to the 
spin-orbit splitting, tend to increase the spin conductiv- 
ity; yet, in the strong disorder regime, the spin splitting is 
notably strengthened as a consequence of the localization 
of states around the impurities. The result is that, in the 
presence of an electric field, the usual torque mechanism 
that drives opposite spins to drift in opposite directions, 
trans ver sally to the electric current, 2 is greatly reinforced 
by the strong local spin splitting. 

Finally, we find that the spin conductivity fluctuations 
increase with the disorder in the extended state regime, 
as already noted in Ref. 9, but decrease above the lo- 
calization transition. Concomitantly, we can follow the 
spin current associated with the wave packet, defined by 
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FIG. 9. (Color online) Static spin Hall conductivity as a func- 
tion of the Fermi energy, for different values of the disorder 
strength. The inset show the variation of a S H at fixed cf — 2; 
in the strong disorder regime it increases with the disorder 
strength. 



jm = mm\m), 

J z (t) = hm J dx [^W T - vivv>j . 

This current vanishes in mean, but undergoes time fluc- 
tuations with amplitudes steadily increasing as a function 
of the disorder intensity, roughly linearly with J s in the 
strong disorder range of parameters (at fixed concentra- 
tion). This behavior is in accordance with the mechanism 
described above of the spin conductivity increase, as it 
is intimately related to the spin current-charge current 
correlations. 

V. DISCUSSION AND CONCLUSION 

We considered a two-dimensional electron gas, with 
spin-orbit coupling and magnetic impurities, to study the 
spin Hall conductivity. We first investigated, using nu- 
merical methods to explore the whole range of parame- 
ters, the dynamical and spectral properties of the quan- 
tum states, as a function of the disorder strength. We 
then studied the corrections to the spin conductivity due 
to quantum interference in the weak disorder regime, us- 
ing the linear response theory. Finally, we computed the 
spin conductivity in the strong disorder regime. 

The evolution of a wave packet shows a rich dynamical 
behavior, with a disorder dependent spreading exponent, 
ranging from almost ballistic to slightly subdiffusive mo- 
tion. This seems to be controlled by the modification in- 
troduced by the spin-orbit interaction, to the simple free 
motion. It is at variance with the usual diffusion observed 




FIG. 10. (Color online) Local and total density of states for 
(left) J s = 1, at the localization transition, and (right) J s = 4, 
in the presence of impurity bands. The (unnormalized) local 
density of states is represented for two arbitrary sites in the 
lattice, occupied by impurities. 



in a random magnetic field, or even in disordered sym- 
plectic systems, and reminiscent to the behavior of quasi 
periodic systems. In addition, the geometry of the wave 
packet, as characterized by the correlation dimension, is 
fractal over the entire range of disorder strengths. In fact, 
we find that for a critical value of the exchange constant 
(at fixed impurity concentration), the system undergoes 
a localization transition. We observe that at this point 
the wave packet width and the correlation dimension are 
compatible with the characteristic diffusion exponents of 
a two-dimensional random walk. Above this point the dy- 
namical inverse participation ratio remains finite at long 
times, indicating the spatial localization of the quantum 
states. This is confirmed by the spectral properties as 
displayed by the local density of states. In the localized 
regime, a strong intermittence arises (arithmetic and ge- 
ometric means do not coincide), characterized be the ap- 
pearance of strong spin correlations between carriers and 
impurities. 

It is natural to think that the interference effects ob- 
served in the phenomenology of the wave packet spread- 
ing, should modify the spin transport. Indeed, the ana- 
lytical computation of the quantum corrections using the 
Kubo linear response theory, shows that these effects con- 
tribute to increase the spin conductivity, and that they 
are proportional to the spin-orbit splitting energy. Yet, 
the effect of the localization on the spin conductivity is 
rather surprising. We found that the spin conductivity 
actually pass a minimum at the transition, and strongly 
increases with disorder in the localized regime. The stick- 
ing of the electronic energy states to the impurities lo- 
cation, and the high spin polarization, contribute to re- 
inforce the spin splitting and enhance the spin density 
fluctuations that, when an electric field is applied, will 
drift in opposite directions to create a strong spin cur- 
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rent. 

It would be interesting to generalize the present one 
particle model in a random potential to take into account 
the electron mediated interactions between magnetic im- 
purities. This would allow us to investigate the interplay 
of Anderson localization and ferromagnetic transitions, 
and their influence on the spin transport, as revealed for 
instance in recent experiments. 56 
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Appendix A: Cooperon 

The total gr-dependent Cooperon can be written in the 
form 



C(q) 











c n \ 
L n ) 





r n 


r n 
u it 








u ti 


u tt 





Wt 











(Al) 



with 



16 



(x 2 + l) 3 ((x 6 + 3x 4 + 2) q 2 l 2 + 2 (x 2 + l) 2 (5x 2 + 4 



u it 



-it 
-it 



-tt 



" u ti 



(O 6 + 3x 4 - 6x 2 + 4) <? 2 / 2 + 4 (x 2 + l) 2 (5x 2 + 4)) ((3x 6 + 9x 4 + 6x 2 + 4) q 2 l 2 + 4 (x 2 + l) 2 (5x 2 + 4)) 

8 e~ 2i Xq 2 l 2 x 2 (> 2 + I) 3 4 + 3x 2 + 6) 
(O 6 + 3x 4 - 6x 2 + 4) q 2 P + 4 O 2 + l) 2 (5x 2 + 4)) ({3x 6 + 9x 4 + 6x 2 + 4) g 2 / 2 + 4 (x 2 + l) 2 (5x 2 + 4)) 

(x 6 + 3x 4 + 2) g 2 / 2 + 2(x 3 +x) 2 



(g 2 / 2 - 4) (q 2 l 2 (3x 2 - 1) - 2 (x 2 + l) 2 (3x 2 + 2)) 
(I0x 4 + 28x 2 - q 2 l 2 (x 4 + 3x 2 + 6) + 26) x 2 + 8 



(g 2 / 2 - 4) (V/ 2 (3x 2 - 1) - 2 (x 2 + l) 2 (3x 2 + 2)) 

where x is the angle between the k and qr and I = tpf is the mean free path length. Because the interaction is 

<tl 

ti 



isotropic on average, Ct\ = after integration over x- The qr-integrated Cooperon is written 



2\2 



(1+X 2 ) 

8ir 2 N F l 2 i 



" lno . f (8+15tt 2 +24s 4 +9tt 6 ) \ 

i0 & ^ (i+x 2 )2(8+i2x 2 ) y 

2 (3x 4 + 8x 2 + 13) 



(4 + 5x 2 ) 



(cr x 



> (J X + (Jy 



> cr^) + (I - cr 2 a 2 ) 



+ 



f T 2 , i \ W f (2+s 2 )(3+5tt 2 +6tt 4 ) \ 

(x 6 + 3x 4 + 2) 



(A2) 



Using the q-integrated Cooperon for the conductivity calculation and nondimensional expressions for the Green 
functions (Eq. (16)) 



6a 



dk 



(G A (k, e F )J z Jk)G R (k, e F )) ® (G R (-k, e F )J x (-k)G A (-k, e F ))C 



= -2nN F (2rf j g- | g (G A (£, <p)j;G R &<p)) ® (G fi (£,^ - tt) J,G^ - tt)) 
5^h = -[ I I (G A (k)J z (k)G R (k)V(0,cj ) )G R (k')) 

JnJkJk' 
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(A3) 
(A4) 



(A5) 
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The expression of the third diagram has the same structure as that of the second diagram. The first quantum 
corrections are written 



S °aH = 



-1 A 



(x 2 +2) 



8 (x 2 + 1) (x 4 + x 2 - 4) 
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The second and third Hikami boxes give an equal contribution 
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